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Composition  Alteration  of  Stratospheric 
Air  Due  to  Sampling  Through  a  Flow  Tube 

1.  INTRODUC  TION 

Siri.'t*  1 H  7  4 ,  ’he  Air  [• '  »r  -r  Goophvsi  s  i.ib>ra’orv  has  been  invoh  ed  ir;  a  re 
s»-  i!'  h  i>r>yra:r-  > ; i*  erned  wi’h  ervoifenh  whale  air*  sampling  a!  The  st  ratnsphore  | 

( .all.iL’ht*  r  an)  Pb-ri  ).  In  :1ns  r  >  r- j  u  ram .  fixed  mer-.il  <  \Tinders  Chree  in  ’.he  ] 

!  • .  » -•*  r»-  » •  r  it  ,  » r  i  ■  i  l!  u  i  a*  l  ms  >,  a  i  r  ii  .sp»-<  iallv  ’  reate.  i  inner  suria.es.  art1  i»v.  tv'-rse  i  j 

! n  a  irjui  1  heliun  !•  i’!i  surrounde  i  bv  a  jfuard  \*>lumc  >!’  liquid  nitrogen.  This  | 

*ri  u  h  be  tir  sar  :>l»-r  r’j  KI\\AS»  :s  she  heart  >f  ’he  balloon  ilijjht  paekaue  sh  twn  j 

a;  l  !  jure  l. 

(/.on  th  >uj!i  *hi-J  program  i  -  rejatjvelv  s*  ra  itjhi !o  rward  in  on-ept,  there  are  ] 

■  .tf,-/  i  s -  iss  trial**  i  will;  ‘he  .  *»lle«  :i  »ji,  r«-(/em-*rati  mi.  and  analvsis  e  nr.  -  ! 

nb-\  >si>h<-r  i.  L«.is  ’■  lxture-  thn*  hav»-  the  n  oten’iai  to  alter  the  jrijinal  ambient 

e  .»  >-i*i  m  »:  ‘lie  -a--  : » i * •  i  l* a"'.  I  he  11  *r»*  impor-tan'  o!'  these  have  been  examine  l  j 

-  '  \ 
r.  '  h*  !.ih  t  )’  >f-v  i  r-  t»l*  • .  -  •  •  ■  <  >  a  !  I  a  l:  he  r  '  ■ '  a  1  “  ).  II  mover,  mil  i  l  n  )  a  .  n  >’  ;  j 

M  \  <*•*•?;* :  m  b.o  r;  ;:vm!  *  >  ‘he  nessihjlitv  of  sample  a  Iterati*)n  i  n  the  relatiw  lv  |j 

l  in.'.  '’]l«  ■  -  'uil'  be!  -  >;•.  u  mi  Ida  >■'  the  lullnn,  a-  sh  twn  in  I'iiiur*-  1.  5  j 

t  Hr-  .•  1  •  ir  :  »■  1 1  *  i .  r  :  *  •  i  :  luu.uv  1  '.'id  )  i- 

t  1 .  ( i  i  I  i.i  j\ i»-  r .  <  .  i  .  .  tt  i  I  h« • ;  ; ,  K .  \  .  t  1  '•  7  m  t  (  rv  >L’eni<  \\  lv>b.‘  Ai  r  San .  f)le  t*  an  i 

I  Pr  *l* !*.**»  ’  >r  st  v  i*  rr  t  >n  i»  >s i :  i  mi  Stu  Ties,  A  1C i  -715 -F)  152 . 

f  >.  ( ia  I  l.iyhf  f,  .....  I  irO»-r-j,  an  I  I  an,  her.  (,.\.  (  1  •  •  M  1  I 

J  St  rat  > spite  i  i .  I  r  i  ♦•  (  ,a  -  (  »:•  ;  >  r  »  >ti  M  u  lies  I  r  i  li  /  i  nt*  I  n  Si  ru  I  rv  >aeni  c  ,  \ 

WVTe-AIr  "I'll  r,  j  \;.  >d-.  \i  (  Ti .  Til  ni  7m7T.‘  ATTTVTfTTTTTT:  ! 


274,000  cutl 

balloon 


38  m 


Figure  1.  Balloon  Flight 
Pai'kagc  for  1RIUAS 


Of  course,  once  the  air  sample  is  removed  fr  nr  'he  ambient  >•  ratos phi •  n- ,  the 
photolvtir  mechanisms  driving  'he  in -situ  steadv -stale  •  hemi-trv  ms  .  !  !«■  re 

sul taut  relaxation  of  reactive  an- 1  semi  -  reactive  . .  ies  during  their  residence  •  i t ■  - •  • 

in  transit  from  the  ambient  atmosphere  to  the  sun  -  pi  i  no  point  has  the  nop-ntta!  to 
moiiit'v  the  initial  relative  compositions  >:  minor  an  i  train  species.  In  a- 1- lit i  at,  the 
eoneent  rat  ions  if  reaetive  species  ami  those  with  sigt'.i'h  eit  probabilities  tor 
hete  rogeneous  removal  tine  to  aiisorpt  i  m  n-.'i-  ti-m  in  the  walls  o!  the  draw  -iti  tube 
could  be  s i uni fi cantl v  reduced  during  their  passage.  1  he  ohirctive  it  this  report 
is  to  assess  the  relative  importanre  if  santpi  tin.’  tube  effects  >n  the  alteration  it 
'he  original  composition  as  sampled. 

There  is  a  strong  resemblance  between  the  operating  elia raete ri sti es  of  the 
sampling  tube  and  those  of  a  tubular,  continuous  flow  chemical  reactor.  I  he  air 
rent  analvsis  makes  use  of  this  analogy  in  its  approach  I  i  developing  a  diagnostic 
model  for  the  sampling  tube.  Kssentiallv,  a  general  -  purpose,  multi -reaction  chett  i 
cal  kinetics  code  KTIKMSF.N)  was  used  to  des.  tube  the  h  miogeneous  g:is  phase 
kinetics,  and,  after  suitable  adaptation,  lu-te rogeneous  wall  effects  as  well.  I  h- 
.  hemical  kinetics  effects  in  the  sampled  air  were  superimposed  on  a  ul  1  v  -  ie\ a -1  one 
laminar  flow  field.  I  lie  development  of  this  approach  and  some  sample  results 
from  the  resultant  model  are  presented  tti  this  report. 


2.  MODEL  DEVELOPMENT 


2.1  Flow  Field  Considerations 

The  Bendwav*  tubing  that  provides  a  constant  flow  of  stratospheric  air  through 
the  sampling  point  (shown  in  Figure  1>  is  about  6  m  in  length  with  an  I.  I  >.  of  7.6  . -r: 
As  shown  in  Figure  2.  at  the  sampling  lorati  >n  three  2.64  -  i:  t).  1).  tubes  i  me  tor 
each  ervo -sampler  volume)  are  located  radiallv,  120°  apart  such  that  thev  ter 
minate  at  the  oircumft  fence  of  a  2.  54  -cm  diameter  circle  .-on.  entr:  with  the  tube. 
Thus,  these  tubes  effectively  sample  the  central  2.  M-cn,  diameter  core  of  the  flow 
in  the  Bendwav  tubing. 


Figure  2.  Schematic  of  the  Crvogemc 
Sampling  Point  Within  the  Sampling  Tube 


Flow  through  the  main  sampling  tube  is  maintained  bv  an  axial  fan  (1  K\\  ,  Tvne 
VAX-2MM  Vaneaxial  blower)  rated  at  40  CFM,  located  downstream,  if  the  sac  pling 
point.  Flight  data  with  the  same  tvpe  of  fan  in  another  program  (see  Sherman'  i 
indicates  that  the  volumetric  displacement  of  this  fan  remains  relativelv  constant 


3.  Sherman,  (lt'82)  private  communication. 


at  the  ra'.i'il  value  until  at)  Hit  30  km,  where  a  reduction  iti  capacity  of  approxi¬ 
mately  am' -tlii  ril  is  indicated.  Thus,  at  rati-. I  capacity,  the  volumetric  flow  rate 

sh  >u 1 , 1  be  about  1.80  ■  10*  cm'* .’see.  which  in  the  80  km  region  drops  to  about 

4  8,  8 

1.  26  -  10  cm  i  see.  I  situ;  these  figures  ami  a  cross  sectional  area  ot  45.  4  cm  . 
'he  average  velocity  in  the  sampling  tubing  is: 

tt  <|  .4  4  IH  cm, sec  (277.5  mi  sec  at  80  km)  . 

4 

Assuming  typical  stratospheric  conditions  (L.S.  Standard  Atmosphere  ): 


15  km: 

2  16.  65 K : 

U0.  85 

torr; 

4.  05 

<io,a 

‘•m 

20  km: 

2  16.  67  K ; 

4  1.  47 

torr; 

1.  8 

<  10 18 

i  •  m 

30  km: 

23  1.  24K ; 

8.  ::o 

torr; 

3.  7 

<  10!7 

rm 

llie  viseositv  of  air  at  these  conditions  (from  Bird  et  al,J  pp«  16-17)is  u  ~  0.015  c  . 
Thus,  the  rorresponding  Kevuolds  numbers  are: 

He  I iSu:  ji 
He  (15  km)  4  110 
He  (20  km)  1827 
He  (30  km)  -  250  . 

Thus,  even  though  the  Kevuolds  number  at  15  km  is  in  the  transitional  flow  regime, 
under  most  sampling  conditions  the  fullv  developed  flow  field  in  the  tube  should  be 
laminar;  that  is.  He  "  2100. 

Taking  the  preceding  into  consideration,  the  conditions  in  the  sampled  central 
core  (-Hj  r  H  j )  will  be  characterized  bv  an  average  velocity,  uKj,  and  an 
average  concentration,  C  .  Averaging  over  the  familiar  parabolic  flow  velo¬ 
city  prol'i  le;  * 


o 


4.  '.S.  Standard  Atmosphere  (1P76)  NOAA,  NASA,  USA  I-',  Washington,  I),  C. 

5,  Bird,  K.B.  ,  Stewart,  W.  K.  .  and  J.ightfoot,  K.  H.  ( I860)  Transport  Phenomena, 

Wilrv  (v  Sons,  New  York. 
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<2l 


1'  if  the  ruiTont  situation,  R  I.  27  cm  aril  R  3.  M  on:,  ana  thus, 

1 


it,,  a  1  -  (  ! .  2  i  )" 

R  i  ’  1 

2(2.  8) 


0.  !*44  tl 


u,,  0.;*44  CJu>  0.1*44  l2'\416>  7  85  rm  sec  (524 

H  1 


■rti  sec  at  30  ktv  ) 


Thus,  the  centra!  core  travels  at  a  mean  velocity  of  735  cm  'sec.  The  mean  core 
oncentration  is  estimate. t  in  Sent  ion  3.3. 

12.2  Axial  Convection  vs  Radial  Diffusion 

In  or.ler  to  assess  the  potential  itn porlance  of  diffusional  losses  to  the  tube  wall, 

the  relative  time  scales  of  axial  eoiivevtion  and  radial  diffusion  must  lie  compared. 

The  general  problem  of  dispersion  in  laminar  tube  fl  nv  has  received  considerable 

6  7  8 

attention  in  the  literature;  the  classic  analvsis  of  lavlor  '  ’  having  been  amplified 

"  10  6 

m  numerous  oeeasions  (for  example,  see  Aids  and  Hunt  ).  Referring  to  1  avion 
the  characteristic  time  for  convection  through  the  sampling  tube  is  given  bv: 


where  1.  is  the  sampling  tube  length  (6  in);  while  that  for  radial  diffusion  is  given  by; 


(  R  ?  3 .  «>“  1 ) 


where  1)  is  the  molecular  diffusivitv  and  K  is  the  tube  radius  (3.8  cm).  A  compari¬ 
son  of  these  time  scales  at  three  altitudes  for  NO  in  air  yields: 

(Due  to  the  large  number  of  references  cited  above,  thev  wall  not  be  listed  here. 

See  References,  page  3f>.  > 


Altitude 
(knt ) 
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2(1 
:i  0 
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■j 

10. 
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'  ii  it 

Is) 

(SI 

0.  7  2 

1.0 

0.  72 

0.  a 

1 .  Oil 

o.  1 

where  the  it:  ilecular  diffusivities  were  estimated  'Ton  Bird  el  al  (Kq,  (16.  U-li, 

5051.  Titus,  at  1  >w  altitudes  7  un  i  7  ..  .  are  of  comparable  magnitude. 

1  .  i fir.  n  h. 

However,  with  it),  reasint’  altim-h  r  ...  r  ,  lue  to  •!•.<•  proportional  increase 

nit.  •  -on  v. 

in  1)  with  tie  ceasing  pressure  (that  is.  Itiv  1  )»»,  )  ,r  purposes  of  estimating  the 
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where  .1  is  'he  zeroth  >r  ier  Bessel  run.  tion  >f  the  first  kin  t,  the  Of  at  e  constants 
i  n 

determined  (iv  the  boundary  ,  ondifions,  A  are  i  'instants  letemine  i  bv  application 

n  '  ■ 

of  orthogonality  of  the  Bessel  functions,  t  is  time,  and  z  -  r  H. 

2.3  Heterogeneous  Wall  Loss  F.stimation  Technique 

In  order  to  use  the  niultireaction  chemical  kinetic  code,  CHfclMKIN  (Kee  et  al1^) 
to  model  the  homogeneous  reaction  behavior  in  the  sampling  tube,  and  also  to  esti¬ 
mate  the  effects  of  wall  termination,  the  heterogeneous  loss  due  to  radial  diffusion 
must  be  estimated  and  transformed  into  an  approximate  pseudo-homogeneous  reac¬ 
tion.  The  development  of  the  resultant  averaging  technique  follows. 


11.  Crank,  .1.  I  1"75)  The  Mathematics  of  I  liiTusion.  Oxford  i'niver.«itv  Press,  London. 

12,  Kee,  H..I.,  Miller,  .1.  A.  ,  and  .Jefferson.  1.11.  <  1  f * B 0 )  CHKMK1N :  A  tieneral- 

Purpose,  Problem,  -Independent .  Transportable,  IOITTHAIX  Chemical  Kinetics 
Code  Package,  SAN  i )  30 -301)3 ,  Sandia  National  Laboratories, 

Livermore,  C  al i fo rn i a . 


At  the  tube  wall  (r  =  H.  or  z  =  1),  a  steady -state  mass  balance  yields; 


-l) 


t)C 


r  -  R 


k  C 


H) 


where  k  is  the  heterogeneous  loss  rate  constant  (cni/sei-l.  Making  use  of  Kq.  (6); 


d  C 

/ 1  \ 

a  c 

r- 

\K) 

'  ITT 

OC 

2  - 1  2 
1) 

' 

\' 

n  - 

A  <y  1 
.  n  n 

iv  .1 

is  the 

first  ordfi 

•  liessel 
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l  n 


pression  into  Kq.  < T > : 


t  1  ;  >  1  ■  2  -  1  > 

D  A  U  "  a  .1 ,  in  “  Id  "it)  cx()i-<v  i> 
n  ■  1  11  1  !1 


n  =■  1 


A  I  uv  1  ’  ~  I )  1  "Cl  exp  t  - (V  (i  . 
n  o  n  n 


Approximating  both  sides  of  12q.  <!>>  In  the  first 
a  better  approximation  as  time  progresses!: 


■) 


I  '  9  -  1  i  ') 

.1  ( cv  -i)  “10 

o  1 

- YT2 — r — rv — 

,1  to-  Id  10 


term 


in  the  series. 


t  which  becomes 


<  101 


The  effective  loss  rate  for  a  particular  species  in  the  central  sampled  core, 
- 1{  |  =  r  ~  Kj  ( K  j  -  1.27  cm),  due  to  steady-state  wall  termination  is  given  by: 


-D  3  C 

TF 


Tt  H, 


2jjU, 


k'  C 


m,  K 


111) 


where  k'  is  a  pseudo-homogeneous  first  order  rate  constant  (sec  '),  and  C  ,,  is 
1  m,  l< 

the  mean  concentration  in  the  sampled  central  core.  liq.  til)  effectively  transforms 
the  diffusive  loss  rate  at  the  boundary  of  the  sampled  central  core  to  an  approxi¬ 
mate,  pseudo-homogeneous  loss  rate  based  upon  the  mean  concentration  in  the  core. 

The  me. in  concentration,  C  ^  is  determined  bv  averaging  the  concentration 
profile  over  r:  * 
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C(i')  r  dr 


'm,  H(  Kj 


l 


r  dr 


(12) 


or,  substituting  Fq.  (6)  for  Civ): 

JO 

K 


C  „=2 

m,  a, 


*11  -  A  J  {a 

C  ,  n  o  n 

/  ^ - 


1  /  9  -1/9 

A  J  (a  ”U  “r)  exp(-a  t)  rdr 

-  '■  n 


« u» 


Solution  of  Fq.  (131,  once  again  using  the  first  term  approximation,  yields: 


C  ,,=2 

nt,  Kj 


12  -  l  > 

.1,  uv.  1.)  K,  i  exp  i  -  a  t> 
11  1  r  ” 


Kj  I  *  1 


- j— t — rp^j — 

a,  U  -  a 


111 


Substitution  of  Fq.  (14>  into  Fq.  (Ill  yields: 


k'  =  a^  . 


l.i  I 


1-or  instantaneous  loss  at  the  wall  (r-  K', 


(.  ..  0  --  -  A  J  (a  1  2DM  “  K>  exp  i  -  a  t) 

r= a  .non  1  n 

n  =  l 
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or,  using  the  first  term  and  first  root  of  J  , 


1/2  -1/2 

(»  j  D  l/“  a=  2.  4 


(  IT' 


Substituting  into  Fq.  (li) 

k’  =  (2.  4'“  D /  K~  i  1111 

which  is  the  approximate  representation  of  the  upper  limit,  pseudo-homogeneous, 
first  order  rate  constant  for  heterogeneous  wait  loss. 

The  rate  constants  representing  wall  loss  can  also  be  interpreted  in  terms  of 
a  fractional  loss  of  the  upper  limit  value  given  by  Fq.  i  181.  i  or  complete  wall 
destruction: 
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\(  wall  flux) 


o  a  c 
a  [■ 


H 


n 
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n=  1 


A  (a 
n  n 


i  n 


- 1  /  •> 

D  R)  J  t« 


-1/2 


l{)  exp  (-a  t)  . 
n 


1 1  to 


Substituting  for  a  j ,  using  Eq.  (17),  and  the  first  term  in  the  summation: 

N  =  DC  2.  4  J  ,(2.  4)/  11 
o  1 


1)  t  1, 25/  II  , 
o 


(20) 


which  represent  the  upper  limit  wall  flux  for  complete  destruction.  If,  f  is  the 
fractional  destruction  at  the  wall,  then  combination  of  Eqs.  (7)  and  (10)  with  the 
definition  of  f  vields: 


f  = 


1/2  -1/2 

kC  .1  (O'  1J  1 '  R) 

o  o  1 

uT  n  257  h 

o 


(21) 


or,  upon  rearrangement: 


k  =  f  D  1.  25/ K  J  (a,  1/2n‘,/2  R)  . 
o  1 


(22) 


Elimination  of  k  between  Eqs.  (10)  and  (22)  results  in: 
1/2 


f  UI/2  1.  25 


1 


R  J  t(0f  1 1  D 


r 77 


R) 


(23) 


Thus,  specification  of  f,  along  with  the  parameters  D  and  K,  defines  via  solution 
of  Eq.  (23).  However,  Eq.  (23)  is  transcendental  in  a  and  thus  requires  an  iterativ 
solution.  An  appropriate  iterative  solution  scheme  using  Newton's  method  is  formu¬ 
lated  in  Appendix  A.  Of  course,  once  a  j  is  determined,  then  k'  follows  directly 

from  Eq.  (15);  namely,  k '  =  a .  A  plot  of  the  heterogeneous  wall  loss  rate  constant 

2  ‘ 

in  dimensionless  form  (k' R  ID)  as  a  function  of  f  is  presented  in  figure  3. 
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Tilt'  corresponding  rate  constant  values  for  specified  f  are: 


f_  iil 

0  0 


0  «.  f  «-  1 


f 


I 


Lq.  i2:i),  k  -  a 
<2.  4) 2  U  U2 


1 


I  igure  i.  Dimensionless 
Heterogeneous  Kate  C  onstant 
as  a  I- unction  of  r ruction  of 
the  Maximum  Wall  l.oss  Kate 


3.  COMPUTATIONAL  COOK 


The  entire  preceding  development  was  undertaken  with  the  objective  of  adapting 
an  existing  code  based  on  CIIKMKIN  (Kee  et  allw),  to  calculate  composition  changes 
along  the  sampling  tube  due  to  chemical  reactions.  The  particular  version  of  this 
code  available  was  intended  primarily  for  the  solution  of  homogeneous  kinetic  prob¬ 
lems  without  transport.  Thus,  the  method  in  Section  2.  3  was  formulated  in  order  to 
transform  the  wall  termination  reactions  into  equivalent  pseudo-homogeneous "  form. 

3.1  CHEMKIN  et  al 

CHLMKIN  is  a  pavkage  of  1  OUTRAN  programs  designed  to  facilitate  general 

chemical  kinetic  modeling  of  multi  reaction,  multispecies  problems.  It  was  developed 

1 2 

at  Sandia  National  Laboratories,  Livermore  by  Kee  and  associates  (Kee  et  al  ). 
Essentially,  it  accepts  the  definition  of  a  chemical  kinetic  mechanism  in  familiar. 
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relatively  "format-free”  chemical  notation,  and  transforms  it  into  1  ORPHAN  code, 
which  can  then  be  manipulated  in  almost  any  manner  the  user  chooses.  1  or  the 
current  applications,  CIIKMKIN  was  used  in  combination  with  a  general  ordinary 
differential  equation  solver,  I.SODI.,  and  appropriate  subroutines  to  transform  the 
FORTRAN  code  from  CHFMK1N  into  the  corresponding  mass  conservation  equations 
for  solution  by  I.SOIJF.  This  particular  combination  was  available  as  SFNsIT,  a 
sensitivity  code  package  developed  by  Kramer  et  al,  *  '  that  features  the  powerful 
new  AIM  (analytically  integrated  Magnus)  modification  of  the  (<i  M  (Green's  function 
method)  for  sensitivity  analysis  of  <  JlJlh  systems. 

3.2  Chemical  Kinetic  Mechanism 

The  chemical  kinetic  mechanism  formulated  fot  use  with  the  (  HI.MKIX  MTNSIT 
code  package  was,  of  course,  essentially  a  basic  stratospheric  c.uemical  kinetic 
model  sans  the  phololvtic  formation  and  destruction  reactions.  The  actual  mechanism 
used  was  extracted  primarily  from  Trie  ■■stratosphere  11*81.  ,  .  ,  *"*  along  with  the 
corresponding  rate  constants  in  appropriate  form  aid  units  for  use  in  the  (  IIKMKIX 
si;  NS  IT  code.  This  mechanism  is  presented  in  I'able  1,  along  with  ttie  corresponding 
rate  parameters. 

In  order  to  model  heterogeneous  wall  'orminstion  reactions,  most  of  the  species 
were  provided  with  a  reaction  of  the  form: 

S  -*  S(  )  . 

The  pseudo-homogeneous  rate  i  on.-Aaiits,  k  ,  for  these  reactions  were  estimated 
according  to  the  method  presented  in  section  2.  3.  Pile  specification  that  appears 
in  the  mechanism  definition  input  next  to  the  1  1  reactions  is  the  negative  of  f,  the 
fraction  of  the  total  wall  destruction  rate  assumed  for  that  species.  I  sing  -f  as 
the  input,  the  subroutine  Kl'KIMI.  'listed  in  Appendix  li >  returns  the  pseudo- 
homogeneous  first  order  rate  constant,  k  ,  for  use  in  the  model. 


13.  Kramer,  M.  \.  ,  Kee,  K.  .J.  ,  and  iObitz,  li.  (1932)  (  IIFMSI.'.N:  A  Computer 

Code  for  Sensitivity  Analysis  of  Llementary  Che  lical  Reaction  ModeHTj 
sA.N'D  82-8^30,  Sandia  National  laboratories,  l.ivcrmorc,  California. 

14.  The  Stratosphere,  11*81,  Theory  and  Measurements  (11*82)  UM(),  NAsA,  I  \\, 

NOAA,  Keport  No.  ITT  NASA*  ^Goddard  Space  1-  light  Center, 

Greenbelt,  Maryland. 


Table  1.  Stratospheric  Kinetic  Model 
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4.  KKSUl.l'S  ANI)  DISCI  SSION 


4.1  Cast’s  and  Data 

Hirer  different  altitudes,  spanning  the  range  tif  tin-  st  im  tosphe  ric  whole-a i  r 
sampling  prog  ram,  .vt  ro  conside  red;  namely,  15,  JO,  a  nil  '50  km.  The.-  initial 
-  oncentratiuns  let’  most  of  (lie  species  Here  obtained  from  tile  Stratosphere  lf.'8l.  .  . 
and  f  ahian  et.  d.  *  ’’  I'iui  sets  of  Ml  and  V  ).,  values  were  used:  averages  of  all  the 
m  id- 1  ditmle  A  It  1 1 .  data,  and  those  of  I  ahian  et  al.  '  ’  Summaries  of  all  values  used 
a  re  presented  in  table  2. 

The  final  f  values  used  for  the  heterogeneous  wall  termination  rent  lions  were 
nmeh  more  diffienlt  to  assign.  Of  course,  values  of  f=  1  are  probab! v  quite 
appropriate  for  all  the  free  radical  species.  However,  tin.-  precise  degree  of  inter- 
adit'ii  of  semi  reactive  am)  stable  species  with  the  Walls  of  tile  sampling  tube  is,  for 
the  most  part,  unknown.  1  ltimatelv,  "f"  values  for  these  species  were  \  a  rted  in 
onler  to  assess  their  effect. 

Cha  raeteri/.ation  of  the  controlling  processes  of  all  tin-  various  heterogeneous 

phenomena  that  can  occur  at  the  walls  of  the  sampling  tube  remains  m  impo"tant 

unresolved  issue.  As  shown  below,  these  heterogeneous  processes  h Me  at  least 

the  potential  to  alter  both  the  rel.  *ive  and  absolute  concentrate  ns  of  the  .Ml  IIV> 

x  x 

species  from  those  in  the  ambient.  The  following  discussion  is  intended  to  highlight 
some  of  the-  potentially  more  important  processes  that  may  occur. 

I  irst  of  all,  c  onditions  at  the  samp) in  Lube  wall  were  assumed  to  be  at  steady- 
state  when  a  sample  is  taken.  This  assumption  is  quite  reasonable  in  \  low  of  the 
fact  that  air  was  continuously  drawn  through  the  sampling  tube,  on.  e  at  altitude, 
for  at  least  a  half  hour  prior  to  the  cryogenic  sample  actuallv  being  taken.  Thus, 
at  the  ambient  temperatures  expected  (see  Section  4.  ID,  the  tuht  wall  is  most 
probably  coated  with  a  layer  of  condensed  water. 

\<  is  known  to  be  rapidly  converted  to  \(  i  b\  re.iction  with  adsorbed  water 
'(Ireenc  and  I’.ist^t.  The  requisite  stoiehiomet r\  is  believed  to  be: 

2M).;  *  ||  O  -*  If.VO  -  l/.XO  .  '2D 


la.  I  ahian,  1’.  ,  Dyle,  .1.  A.  ,  and  Wells,  It.  .f,  <  I  b!>2i  Diurnal  vai’ia.if>ns  of  minor 
constituents  in  the  stratosphere  modeled  as  a  function  of  latitude  and 
season,  .1.  (leophvs.  lies.  f!7:4hfU. 

—  —  —  v  VM 

1  .  (iffM-nu,  S.A.  ,  and  Pust,  li.  (  Determination  <>f  nitrogen  dioxide  hv 

solid  chromatography ,  Ana  1.  (  hem.  S0:l0:*o. 
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Table  2.  Initial  Species  Concentration 


Species 

Concentration 

1  mol  fraction) 

Source 

1  a  km 

20  km 

20  km 

No 

1  e .  d ) 
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I  nde  conditions  where  nitrous  arid  is  unstable  (that  is,  at  ruuni  temperature  and 
above),  it  would  decompose  according  to  tin*  stoichiometry: 
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efficiently  trapped  by  condensed  water  (most  probably  as  mixed  clathrate  hydrates! 
at  relatively  high  temperatures,  exactly  in  the  range  expected  for  the  sampling  tube 
walls  (that  is,  21G  to  230  Kb  Rather  than  being  true  thermodynamic  condensation, 
this  phenomenon  is  a  complex  kinetic  entrapment  process  favored  by  high  fluxes  of 
both  the  species  and  water  vapor.  In  the  case  of  the  sampling  tube,  however,  even 
though  the  temperatures  at  altitude  are  in  the  requisite  range,  the  fluxes  seem  to  be 

1  it  20  -2  -  1 

too  low  for  this  effect  to  be  important;  that  is,  10  -  10  m  sec  for  botli  \<  > 

H  IV  -  2  -  1 

and  >  in  the  experiments  cited,  vs  ~  3  ■  10  and  ~  2  •  10  m  sec  ,  re¬ 
spectively,  at  the  conditions  in  the  sampling  tube  case.  Based  mi  these  considera¬ 
tions,  it  is  judged  that  \c>  is  most  probably  not  removed  to  any  significant  extent 
at  the  sampling  tube  walls;  that  is,  f  0,  would  he  a  reasonable  assumption 

Although  the  preceding  discussion  is  intended  io  focus  on  the  most  like!',  fimt 
order  effects  expected  to  occur  at  the  tube  wall,  it  is,  perhaps,  more  re.,|i-ti.  to 
think  of  the  behavior  of  the  tube  at  high  altitude-,  where  radial  diffusivitie -  .re  f.rge, 
as  approximating  that  of  a  chromatographic  column.  The  intern,  tens  m  ..mg  the 
various  species  and  the  condensed  phase  on  tile  tube  wall  can  In-  quite  complex  ,uh 
nonlinear,  and  thus,  almost  impossible  to  estimate  .  prion.  (  ui  reiith,  the  .ini', 
manner  in  which  to  obtain  realistic  "f"  values  for  the  various  spe.  ies  under  -nmplmg 
conditions  ami  concentrations  is  experimentally  . 

4.2  Model  Results 

■4.2.1  c;  i:\KKAl. 


Some  model  results  for  various  cases  of  heterogeneous  wall  loss  arc  ore  -amt 
m  Table  3  |(at,  (hi,  and  (e)  for  la,  20.  .ml  30  km.  respe.  tlx  eh  |,  in  term.-  ..I  fr.. 
'ions  of  initial  ambient  concent  rations  at  'lie  final  ervogenic  sae  pling  muni .  11m 

conditions  for  the  various  cases  are  included  in  the  legend  for  I  aide  :■;.  1  h<-  -an 


results  are  presente.l  in  Table  4  |<a>,  (h),  and  (. 
1  he  cases  in  both  Tables  3  and  4  are  the  same. 
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1  able  3.  Rati  is  if  Species  ( '  nnvnt  lat  t  ms  a*  the  Satruling  Point 
(at  la  km,  lid  20  km,  and  let  30  km 
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<  ase  l.  f  0  l'or  species 

l  a  se  1 1 .  f  1  fo  r  a  1 1  speeierf  exeept  1 3  ),  (  (  )  ,  (  l.j,  \  a  ml  1 1  i  f  0 1 . 

l  ase  III.  f  0.  a  for  .11  species  except  (it,  t  (  )  ,  (  t  \  f  .ml  ||i;  .f  m. 

I  ase  1  \  .  i  -  1  for  all  species  except  \<  ),  ft),  (  (  \.)t  <  1  ,  \  md  II  <f  01, 

I’ he  fol  lowing  fou  r  eases  a  re  similar  to  the  p  re.  .ding  except  'hat  V  >  ..ml  V 


are  taken  from  1  abian  et  al. 

(  ase  \  .  equivalent  to  Case  I. 

Case  \  i.  equivalent  to  C,.se  II. 


!  .qmva  lent  f.  i  (  ,.  -e  III. 
I  .uui  va lent  to  (  ase  1  \  . 


Table  3.  Ratios  of  Species  Concentrations  at  me  Sampling  Point  to  U,e  Ambient, 
(a)  15  km,  (b)  20  km,  and  (c)  30  km  (Contd) 


Species 

Case 

(a)  15  km 

I 

11 

III 

1\ 

V 

V 1 

\  11 

V  111 

NO 

0.  995 

0.  737 

0.  9  1  9 

0.  996 

0.  996 

0.  7  36 

0.  9  l  6 

0.  996 

no2 

1.  001 

0.  836 

0.  931 

0.  76  1 

1.  004 

0.  762 

0.  934 

0.  763 

no3 

23.  775 

12.  535 

21.315 

15.  6  3 5 

16.  152 

16.  000 

21.  7  90 

15.  998 

HN°3 

1.  004 

0.  845 

0.  937 

0.  773 

1.  000 

0.  77  1 

0.  934 

0.  77  1 

UNO, 

4 

1.  000 

o 

£ 

oc 

0.  936 

0.  778 

0.  999 

0.  778 

0.  9  35 

0.  777 

N2°5 

1.  000 

0.  854 

0.  938 

0.  7  85 

1 .  000 

0.  7M  5 

0.  938 

0.  7  85 

CO 

1.  000 

1.  000 

1.  000 

1. 000 

1 . 000 

1.  000 

1.  000 

i .  000 

l'(>2 

1.  000 

1.  000 

1. 000 

1. 000 

1 . 000 

1 .  000 

1.  000 

1.  000 

(  ) 

'0. 

'0. 

-o. 

-0. 

'0. 

'0. 

•o. 

'0. 

°2 

1.  000 

1 .  000 

1 .  000 

1  .  000 

1 .  000 

1 . 000 

1.  000 

1.  000 

°3 

t.  000 

0.  837 

0.  9  -l  1 

0.  76  1 

1 .  000 

0.  76  1 

0.  93  1 

0.  7  6  1 

H(  i 

0.  896 

0.  74  2 

0.  8  1  8 

0.  637 

0.  "64 

0.  t>oM 

0,  Mo 3 

0 .  ti  3  M 

H<>, 

1.  252 

0.  955 

1 .  15  1 

0.  "0, 

1.216 

0.  95  9 

1.21, 

0.  9  5" 

U2°2 

1.  000 

0.  825 

0.  92."> 

0.  7  4  *5 

1  .  000 

0.  746 

0.  892 

U.  7  )■> 

H2° 

1.  000 

0.  7  97 

0.  9]  2 

0.  .Tin 

1 . 000 

0.  7  06 

0.  57" 

0.  7  06 

N2 

1.  000 

1 .  000 

1.  000 

1 .  000 

1.  000 

1 .  000 

1 . 000 

: .  ooo 

CIONO, 

1.  000 

0.  852 

0.  937 

0.  7  M3 

1.  000 

0.  783 

0.  93, 

0.  7  M3 

lit'l 

1.  000 

0.  828 

0.  "27 

0.  750 

1 . 000 

0.  750 

0.  9J7 

0.  7  50 

Cl 

12.  280 

10.  130 

1  1.  350 

1  1 .  960 

0.  3  24 

0.  24  5 

0.  305 

0.  i  2 2 

(TO 

0.  99 1 

0.  831 

0.  923 

0.  7  "3 

0.  ""9 

0 .  7  6  3 

0.  930 

0.  508 

II 

0.  103 

0.  084 

0.  0"4 

0.  07  3 

0.111 

0.  07  7 

0.  099 

0.  0  7  7 

IK  XT 

1.  000 

0.  839 

0.  "32 

0.  765 

1 .  000 

0.  M3 9 

0.  9  3  2 

0.  7  o5 

H, 

1.  000 

1 . 000 

1 .  000 

1 .  000 

1.  000 

1. 000 

1.  000 

1 . 000 
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111 
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VI 
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\  III 

\<  i 

0.  99  1 

0.  5  1  7 

0.  834 

0.  008 

0.  004 

0.  5  1  7 

0.  834 

0.  993 

vij 

1 . 005 

0.  554 

0.  850 

0.  555 

1. 002 

0.  553 

0.  857 

0.  5  54 

NO, 

> 

2.  150 

1 . 056 

1. 77  85 

1. 04  1 

1.  2071 

7  0 .  b !  1 7 

1.  109 

0.  6  96 

H\t  i.. 

a 

1.  000 

0.  57  0 

0.  862 

0.  570 

1. 000 

0.  57  0 

0.  862 

0.  570 

ii\<  > 

1 . 000 

0.  580 

0.  866 

0.  580 

1. 000 

0.  580 

0.  866 

0.  580 

N2°5 

1. 000 

0.  502 

0.  87  1 

0.  502 

1. 000 

0.  592 

0.  866 

0.  592 

(.'<  > 

1 .  000 

1.  000 

1. 000 

1. 000 

1. 000 

1.  000 

1.  000 

1.  000 

rn2 

1.  000 

1.  000 

1.  000 

1.  000 

1.  000 

1.  000 

1.  000 

1.  000 

(  ) 

'0. 

0. 

'0. 

'0. 

'0. 

'0. 

-o. 

'0. 

U2 

1.  000 

1.  000 

1.  000 

1.  000 

1. 000 

1. 000 

1.  000 

1.  000 

( 

.1 

1. 000 

0.  554 

0.  85b 

0.  554 

1.  000 

0.  554 

0.  856 

0.  554 

IIO 

1.  124 

0.  5  1  2 

0.  0  1  2 

0.  540 

0.  965 

0.  446 

0.  776 

0.  447 

HO., 

0.  068 

0.  5  15 

0.  819 

0.  5  1 0 

1. 003 

0.  530 

0.  848 

0.  530 

U2°2 

1. 000 

530 

0.  846 

0.  530 

1.  000 

0.  5  30 

0.  846 

0.  530 

11,0 

1. 000 

0.  47  1 

0.  820 

0.  47  1 

1.  000 

0.  47  1 

0.  820 

0.471 

X2 

1.  000 

1.  000 

1.  000 

1. 000 

1. 000 

1.  000 

1.  000 

1.  000 

OlONO, 

1.  000 
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0.  869 

0.  588 

1. 000 

0.  588 

0.  86  9 

0.  588 

IK  i 

1. 000 

530 

0.  848 

0.  536 

1.  000 

0.  5  36 

0.  84  3 

0.  536 

ci 

14.  007 

7.  845 

12.  582 

14.  430 

0.  622 

0.  328 

0.  526 

0.  604 

CIO 

0.  0.06 

0.  557 

0.  854 

0.  555 

1. 000 

0.  558 

0.  857 

0.  558 

n 

0.  007 

0.  003 

0.  006 

0.  004 

0.  006 

0.  003 

0.  005 

0  003 

IKK  1 

1. 000 

0.  560 

0.  858 

0.  560 

1.  000 

0.  56  0 

0.  858 

0.  560 

"2 

1. 000 

1. 000 

1. 000 

1. 000 

1.  000 

1 .  000 

1.  000 

1.  000 
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Table  4.  NO'NOx  for  the  Eight  Computational  Cases  in  Table  3, 
(a)  15  km,  (h)  20  km,  and  (c)  30  km 


(a)  15  km 

Ratio 

Ambient 

1 

11 
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IV 

NO  No  'a> 

0.  4  17 

0.  300 

0.  274 

0.  298 

0.  359 

x 

N<  >  (N(  >-  NO./ 

0.  4  17 
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0.  386 

0.  4  14 
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N( )  N( 
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\  1 
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\  III 
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Table  4.  NO/NO*  lor  the  Eight  Computational  Cases  in  Table  3 
la)  15  km,  (b)  20  km,  and  <c)  30  km  (Contd) 


( e )  30  km 


Ratio 

Ambient 

1 

11 

111 

IV 

NO/NO  u° 

X 

0.  248 

0.  244 

0.  168 

0.  195 

0.  967 

N()/(.\’()4  no.,) 

0.  345 

0.  339 

0.  260 

0.  283 

0.  981 

NO/. MO 

0.  528 

0.  513 

0.  352 

0.  395 

0.  525 

Ambient 

V 

\  1 

V!1 

\  111 

NO/NO  (a) 

X 

0.  080 

0.  07  9 

0.  042 

0.  068 

0.  880 

NO/1NO+  NO.,) 

0.  167 

0.  164 

0.  107 

0.  149 

0.  95  2 

no/.no2 

0.  200 

0.  195 

0.  120 

0.  175 

1  9.  920 

<a,.NO/NOx  is  taken  as  NO /(NO  •  NO.,  •  NO.j  ■  HNO.^  •  H\04  -  2\.,0  >. 

Cases  1  and  \  in  Tables  3  and  4  show  the  effects  of  the  homogeneous  gas  phase 
reactions  only;  that  is,  no  heterogeneous  wall  loss  for  any  of  the  species  if-  O'. 

The  behavior  of  species  in  the  sampling  tube  in  these  tivo  eases  is,  for  the  most 
part,  predictable.  Due  to  the  relatively  short  average  residence  time  of  the  sampled 
core  (0.76  sec  at  15  and  20  km,  and  I.  14  sec  at  30  km1,  only  those  reactions  with 
time  constants  comparable  to  or  less  than  the  residence  time  have  any  appre¬ 
ciable  effect.  Consequently,  the  most  dramatic  change  in  the  ambient  composition 
is  the  rapid  decrease  of  ()  atom  concentration  upon  sampling,  due  to  reaction  25 
(see  Table  l).  The  half-life  of  O  in  the  sampling  tube  due  to  this  reaction  varies 
from  1.7  3  10  ^  sec  at  15  km  to  0.  024  see  at  30  km.  The  time  scales  of  most  of 
the  other  reactions  listed  in  Table  1  are  significantly  greater. 

At  15  km,  where  diffusional  losses  to  the  sampling  tube  wall  are  relativelv 
unimportant,  the  ambient  composition  is  not  seriously  altered  under  almost  any  set 
of  assumptions  regarding  wall  loss  "f"  values  (see  Table  3a),  especially  with  re¬ 
spect  to  those  species  related  to  the  oxides  of  nitrogen.  This  is  also  reflected  in 
the  fact  that  the  NO/NC>x  ratios  in  Table  4a  remain  relativelv  constant  for  the  eight 
computational  cases  considered.  At  20  km  (Tables  3b  and  4b),  the  potential  effect 
of  diffusional  losses  to  the  tube  wall  becomes  evident,  and  at  30  km  (Tables  3c  and 
4c)  the  effect  can  be  quite  severe.  However,  the  relatively  small  sensitivity  of  the 
absolute  extent  of  wall  loss  to  species  identity  shows  that  the  variation  in  diffusivitv 
among  species  is  clearly  a  second  order  effect.  Also,  as  shown  bv  C  ases  1  and  \ 

(f=  0)  in  Tables  3a,  3b,  3c,  the  homogeneous  gas  phase  chemistry  in  the  sampling 
tube  is  generally  incapable  of  appreciably  altering  the  sampled  composition  from 
that  in  the  ambient;  the  residence  tinies  are  simply  too  short. 
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Of  the  species  examined,  VO  and  Cl  exhibit  large  increases  in  concentration 
as  a  result  of  the  homogeneous  gas  phase  chemistry.  NO.^  is  observed  to  increase 
appreciably  in  Tables  3a  and  3b  due  to  production  from  NOg  oxidation  by  ozone 
(that  is,  reaction  6  in  Table  1).  However,  N'O^  is  also  rapidly  consumed  by  NO 
(that  is,  reaction  7  in  Table  11.  Thus  the  ratio. 


k  |0  ]  (NO  ] 
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controls  the  behavior  of  NO.(  in  the  sampling  tube  due  to  homogeneous  gas  phase 
chemist  y.  In  addition,  as  shown  in  Table  3c,  the  \t>3  concentration  can  be 
severely  affected  bv  heterogeneous  wall  loss,  to  the  point  where  production  bv  V  > 
oxidation  is  completely  overwhelmed,  such  that  VO ^  decreases  Rionotonicall.v  with 
progress  through  the  flow  tube.  In  any  case,  the  No^  behavior  should  not  nave  a  in¬ 
appreciable  effect  on  tlie  total  VO  anahses,  due  to  its  relatively  low  expected 
concentration. 

The  behavior  of  Cl  is  somewhat  analogous  to  that  of  NO  .  As  can  be  seen  from 
Table  3a,  in  Case  1  the  Cl  concentration  increases  by  about  a  factor  of  twelve, 
while  for  Case  Y  it  actually  decreases  by  approximate!-,  r.-.o- thirds,  i  ron;  the 
legend  of  Table  3,  the  only  difference  between  Cases  1  and  11  is  the  absolute  values 
of  the  ambient  concentrations  assumed  for  NO  and  N<  >.,.  obviously,  in  Case  1,  Cl 
increases  due  to  the  reaction  of  do  with  No  i  that  is.  reaction  -1  in  Table  !>;  while 
in  Case  V,  Cl  decreases  due  to  the  reaction  of  cl  a.  ith  <  > ,  uhat  is,  reaction  20  in 

a 

Table  0.  Thus,  the  primary  result  of  whether  Cl  increases  or  decreases  during 
passage  in  the  sampling  table  depends  on  the  ratio: 


k4  [NO]  | CIO] 
k20  (Cl]  [03|  • 


At  20  km  the  effect  is  quulitativ  lv  the  same  as  for  15  km.  At  30  km,  however.  Cl 
decreases  monotonieallv  from  the  ambient  for  both  Cases  1  and  \  .  The  primary 
difference  at  30  km  is  that  the  heterogeneous  wall  loss  becomes  significant  enough 
to  overcome  Cl  production  via  reaction  4,  for  the  No  concentrations  assumed. 

This  sensitive  behavior  of  Cl  suggests  an  interesting  experimental  approach 
for  either  NO  measurements  or  the  assessment  ol'  heterogeneous  wall  loss  rates. 
Accurate  values  of  the  rate  constants  and  i  1  measurements  as  a  function  of  distance 
in  a  sampling  tube  >  using  resonance  fluorescence,  for  example'  could  be  used  to 
determine  either,  knowing  something  about  the  other. 
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4.2.2  NO  / HNO 
x  x 

As  concluded  in  the  preceding  section,  and  as  is  evident  upon  comparison  of 

the  three  different  NO/NO  ratios  for  Cases  I  and  V  with  their  respe  live  ambient 

x 

values  (that  is.  the  ratios  determined  from  the  ambient  concentrations  assumed  at 

the  mouth  of  the  sampling  tube)  in  Tables  4a.  4b.  and  4c,  toe  homogeneous  gas 

phase  chemistry  in  the  sampling  tube  is  incapable  of  appreciably  altering  the 

NO  /HNO  sample  composition.  Thus,  the  onlv  manner  in  which  the  sampling  tube 
xx 

can  affect  the  sampled  ambient  composition  is  via  heterogeneous  interae'  j  ons  with 

the  tube  walL.  If  this  is  the  case,  then  the  m  >st  severe  effects  on  sample  i  c  impost  ■ 

tion  should  occur  at  the  highest  altitudes  where  the  radial  diffusivities  are  greatest. 

Examination  of  the  30  km  results  presented  in  Table  4c.  reveals  that  of  the  cases 

studied,  only  IV  and  \  III  result  in  appreciable  alteration  of  the  three  NO  NO  ratios 

from  the  ambient.  Reference  to  the  legend  in  Table  4  reveals  that  in  these  cases 

NO  is  assumed  not  to  be  lost  at  the  wall  (that  is,  f  =  0>,  while  all  the  other  NO  /HNO 

x  x 

species  are  assigned  f=  1.  As  discussed  in  Section  4.  1,  this  scenario  is  not  totally 

unreasonable  in  view  of  what  is  known  about  the  behavior  of  these  species,  and 

would  explain  significant  relative  enhancement  of  NO  over  all  the  other  NO  ''UNO 

species,  as  compared  to  ambient  values.  However,  in  order  for  this  behavior  to  be 

predominant  in  the  sampling  tube,  the  magnitude  of  its  effect  must  increase  with 

altitude,  and  this  is  at  variance  with  the  data. 

In  general,  the  AFGL  stratospheric  NO  and  NO.,  mixing  ratios,  determine.!  bv 
1 1> 

Gallagher  et  al.  seem  to  be  slightly  lower  than  some  other  experimental  observa¬ 
tions  at  high  altitudes,  and  slightlv  greater  at  low  altitudes.  The  fact  that  the 
absolute  values  are  of  the  same  order  of  magnitude  as  those  determined  bv  other 
researchers  with  different  experimental  techniques,  indicates  that  severe  wall 
losses  of  NO  and  NO.,  at  high  altitudes,  as  exemplified  bv  Case.'  11  and  VI  in 
Tables  3a,  3b,  and  3c  (that  is,  for  f  1),  are  improbable.  Also,  the  severe  dis¬ 
crimination  against  NO  in  comparison  to  NO,  as  exemplified  bv  Cases  IV  and  \  III 
in  Table  3c  (30  km,  f  1)  for  NO  /HNO  and  f  0  for  NO),  seems  to  be  just  as  im¬ 
probable.  However,  a  comparison  of  "averaged"  30  km  v  alues  reveals  that  the 
AFGL,  results  are  roughly  one-third  of  some  of  the  other  experimental  results. 

Thus,  if  heterogeneous  wall  loss  in  the  sampling  tube  was  the  sole  course  of  this 
discrepancy,  then  an  overall  average  value  of  f -  0.  5  wall  loss  (that  is.  Cases  III 
and  VII  in  Table  3c)  would  explain  the  results  relatively  well.  If  this  assumption 
applied  at  20  km,  the  sampled  NO  and  NO^  values  would  be  83  and  86  percent  re¬ 
spectively,  of  the  ambient  values  (that  is.  Cases  III  and  VII  in  Table  3c),  which 
would  place  them  in  approximate  agreement  with  other  experimental 


19.  Gallagher,  C.C.,  Forsberg,  C.A.,  and  Pieri,  R.V.  ( 1983)  Stratospheric 
^O,  CFgCIg,  and  CFCl^  composition  studies  utilizing  in  situ  cryogenic 

whole  air  sampling  methods,  .J.  Geophys.  Res.  88:3798. 
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observations  (assuming,  of  course,  that  the  other  observations  have  no  inherent 
experimental  bias  with  respect  to  sampled  vs  ambient  values!.  The  actual  case, 
however,  is  that  the  AFG1.  20  km  data  are  approximately  a  factor  of  two  to  three 
times  greater  than  other  experimental  values. 

Thus,  there  seems  to  be  no  consistent  pattern  in  the  Al'Gl.  NO  and  NO,,  results 
that  unequivocally  implicate  heterogeneous  wall  interactions  as  the  primary  factor 
contributing  to  experimental  bias  introduced  bv  the  sampling  tube  of  the  cryogenic 
whole  air  sampler.  However,  there  remain  some  more  speculative  possibilities 
that  seem  to  be  consistent  with  both  significant  heterogeneous  effects  and  the  ob¬ 
served  data.  For  example,  there  is  the  possibility  that  the  lower  altitude  data, 
which  are  relatively  free  of  wall  effects,  are,  therefore,  more  accurate  representa¬ 
tions  of  ambient  values.  In  this  case,  the  higher  altitude  data  might  be  significantly 
affected  bv  wall  loss,  and  thus  be  biased  towards  lower  values.  Gallagher  et  al“^ 
present  some  results  of  a  more  recent  two-dimensional  stratospheric  model  pre¬ 
diction  (Sze  et  al“  1 )  that  seem  to  suggest  this  behavior;  that  is,  the  model  predicts 
higher  NO  and  NO,,  values  than  the  AFG1.  data  at  high  altitudes.  (However,  it  also 
predicts  lower  values  than  the  AFG1-.  data  at  low  altitudes.  )  Of  course,  if  this 
scenario  had  any  validity,  this  would  implv  that  other  experimental  measurements 
are,  in  effect,  too  low.  However,  this  possibility  is  not  totally  inconceivable, 
given  possible  experimental  biases  inherent  in  the  other  techniques,  and  natural 
local  and  temporal  variations  of  NO  and  NO^  mixing  ratios;  either  or  both  of  which 
can  account  for  the  lower  values  observed. 

Another  possibility  that  is.  in  a  sense,  the  inverse  of  the  preceding  one,  and 
which  also  seems  consistent  with  the  known  data,  focuses  on  the  operational  charac¬ 
teristics  of  the  balloon-borne  orvogenic  whole  air  sampler.  Kssentiallv,  all  the 
ervogenie  whole  air  samples  were  taken  on  the  descent  leg  of  the  balloon  flight  in 
order  to  avoid  contamination  from  outgassing  of  the  balloon  and  gondola.  However, 
this  procedure  exposes  the  sampling  tube  wall,  at  relatively  low  temperatures,  to 
t lie  highest  ambient  concentrations  of  NO  and  NO,,  at  high  altitude  prior  to  collection 
of  the  samples.  Kxposure  at  high  altitudes  occurs  for  enough  time  at  conditions  of 
high  radial  flux,  due  both  to  high  radial  diffusivities  and  high  ambient  mixing  ratios, 

lo  accumulate  significant  amounts  of  NO  /HNO  ,  and  even  NO,  if  the  actual  capture 

xx 

coefficients  are  appreciable.  The  descent  of  the  balloon  to  lower  altitudes  creates 
an  unsteady-state  situation  with  respect  to  wall  accumulation  due  to  decreased 


20.  Gallagher,  G.  C.  ,  Forsberg,  C.  A.  ,  Fieri,  K.  V.  ,  and  Kaueher,  G.A.  (1983a) 

Oxides  of  Nitrogen  Content  of  Whole  Air  Samples  Obtained  at  Altitudes 
From  12  to  30  km,  submitted  for  publication. 

21.  Sze.  N.  1).  ,  Kn,  M.K.  VV.  ,  Livshits,  M.  ,  Wang,  W  .  C.  .  and  Hvan,  F.  B.  (1982) 

A  Heseareh  Program  for  Atmospheric  Chemistry,  Kadiation  ami  Dynamics, 
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I  radial  fluxes  of  NCWHNO^  from  tlie  gas  phase  as  a  result  of  both  lower  mixing 

!  ratios  and  lower  radial  diffusi vities.  For  example,  NO  and  NOg  radial  wall  fluxes 

near  the  sampling  tube  mouth,  estimated  from  Fq.  (19)  (that  is,  for  f-  1)  using 

1 5 

ambient  mixing  ratios  from  Fabian  et  al,  at  the  three  altitudes  considered  here  are 


Altit  ode 

N(NO) 

mno2> 

(km ) 

,  ~2  -I, 

(cm  sec  ) 

,  -2 
(rm  sec 

20 

1.  2  X  10!‘ 

6.  2  <  10'1 

20 

1. 2  <  1  O' 1 

.1.4  X  10H 

15 

5. 2  <  10’ 

4.  8  X  10‘ 

As  the  wall  flux  of  these  species  Iron  the  gas  phase  increases  during  descent,  sort  e 

if  t he  accumulated  NO  UNO  vx  i  1 1  tend  to  evunorate  due  to  the  prevailing  unsteady  - 
x  x 

state  situation  ereated,  therein’  producing  a  net  flux  of  species  into  the  sap  nled  gas; 
that  is,  tile  tube  walls  effectively  outgas.  '1  his  effect  would  tend  to  itvrease  the  NO 
and  NO.,  values  over  those  in  the  an  bierit  for  all  lower  altitudes  during  the  course 
of  a  particular  flight  until  the  wall  accumulation  has  been  significantly  lissinate  I. 

This  hypothesis  qualitatively  explains  son  <■  o'  the  salient  characteristics  of 
die  AFfiL  data.  First  of  all,  due  to  die  expected  low  tube  ■  all  temperatures  (see 
Section  4.2),  the  net  effusion  fr-ot;  the  walls  at  lower  altitudes  would  be  relatively 
slow.  This  would  tend  to  account  for  'he  small  absolute  differences  between  the 
NO  and  NO,  mixing  ratios  determined  bv  die  \FC.1.  r:  (•asurements  and  tho-e  fron 
other  techni que s .  An  oxar  i nation  of  tin  Al  til.  data  reveals  that  although  the  N(  ) 
values  exhibit  a  drop -off  w  ith  decreasing  altitude,  the  -otal  NO  -  NO  results  ire 

>  fj 

remarkably  constant  as  a  function  ol  altitude  (Gallagher  et  al"  I.  This  n  av  in  pjv 
a  conversion  of  NO  to  NO,  in  the  condensed  phase  on  the  tube  wall  aceording  to  the 
well-known  stoichiometry  (for  example,  see  (  hi  lion""  h 

NO  NO.,  11,0  -•  J1INO.,  .  (J  o 

The  nitrous  ai  id,  which  i  -•  unstable  in  the  gas  phase,  would  then  lecon  pose  upon 
desorption  it  a  lower  altitude  ae.ording  to  Fq.  (2a  I; 

till  NO.,  -  HNO  JNO  11,0  .  (25) 


22.  Chilton,  4.11.  (F'fi't)  Strong  Water.  Nitre  Acid;  Sour.es,  Methods  of 
Manufaeture,  and  (5ses  ,  XT77  TVo s s ,  T'anTTri  Tl'^ ,  TiT;* s s .  Tn  j s r* n . 


The  net  effect  of  liqs.  (29)  and  (25)  is: 


HNO.,  -  NO.,  -*  HNO.^  *  NO  <301 

nr,  reduction  of  NO.,  to  NO  on  a  one -for-one  basis,  t  hus,  although  NO  would 
decrease,  the  total  NO  •  NO.,  would  remain  relatively  constant,  since  NO.,  is 
produced  at  the  expense  of  NO. 

As  a  rough,  order-of -magnitude  estimate  of  such  a  process,  assume  that  bv 

whatever  mechanism,  the  tube  walls  become  saturated  with  NO.,  at  20  km  due  to 

pre-sampling  exposure.  Upon  descent  to  20  km,  the  unsteady -state  situation  created 

thereby,  would  produce  a  net  flux  into  the  gas  phase  of  approximately 
n  B  h  -•)  -I 

(6.  2  '<  10  -  3,  4  *  10  I  -  5.  86  ■  10  cm  “  sec  .  For  a  tube  wall  surface  area  ot 

4  2 

600  cm  \  (7  i.6  cm)  -  1.43  •  10  cm  ,  the  net  molecular  How  into  the  sampled  gas 
would  be  8.4  ■  lo'^  molecules  si  c.  Using  a  mixing  ratio  of  3  •  10  at  20  sr 

(Fabian  ot  al'*).  and  the  number  density  and  volumetric  flowrate  in  the  sampling 
tube,  the  total  inflow  of  NO.,  under  these  conditions  would  be; 

(3  <  10 ' 1  °)  (1.  8  •  1018  cm  '  S  (1.  80  -.  K)1  cn  *  sed-  1  •  10  M  molecule.-;  Se,  .  (  )! 

course,  the  estimated  outgassing  rate  of!'..  4  •  In'’  molecules  -m.  will  tend  to 
decrease  with  time,  its  instantaneous  value  determined  m  a  .  m  plex  ta-diion  bv 
in-flight  operating  faetors  and  tube  wall  exposure  hislorc,  and  it  will  not  be  .  on  ■ 
stant  over  the  entire  sampling  tube  length  as  assumed,  so  this  i<  mlv  a  •,  erv  at.nroxi 
male  value.  Nevertheless,  the  obvious  impheatiou  of  this  estiu  ate  is  that  the  out 
gassing  rate  can  be  of  the  same  order-of- magnitude  as  the  sampled  ambient  rate, 
thereby  lending  some  credence  to  this  hypothesis. 

This  "outgassing’  hypothesis  mav  also  be  consistent  with  the  A  Id ,  1 ,  rvo- 
sampler  data  on  nitrous  oxide  (N.,0)  and  Fluoroearbons  11  and  12  (Gallagher  et  ai'  i. 
These  species  are  known  to  be  stable  and  their  mixing  ratios  decrease  with  aitrude 
in  (he  stratosphere.  Thus,  the  expected  effect  of  the  tube  wall,  consistent  with  the 
preceding,  would  be  qualitatively  quite  different.  At  high  altitudes  the  total  wall 
flux  would  tend  to  decrease  with  altitude,  primarily  due  to  lower  at;  Incut  ■  :\mg 
ratios.  Therefore,  the  sampling  tube  wall  would  experience  increasing  wall  :'!u\e- 
during  descent,  which  effectively  inhibit  outgassing.  Thus,  for  these  inert  see,  je-, 
tile  only  possible  role  of  die  tube  wall  is  as  an  adsorber  or  sink,  which  w  mi  l  tend 
•o  decrease  sampled  compositions  from  those  in  the  ambient  at  high  altitudes  where 
radial  di f fusivities  are  large.  Although  this  trend  is  evident  in  the  data  presented 
bv  Gallagher  et  al '  for  N.,0,  and  I*  luoroearbons  11  and  12,  when  eomnare  l  w  ith 
some  model  profiles  (that  is,  good  agreement  al  lower  altitudes  with  increasing 
short  fall  at  higher  altitudes),  it  is  bv  no  means  unequivocal  since  the  data  also  -how 
extremely  good  agreement  with  other  model  predictions  (see  Gallagher  et  al'  ). 
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Of  course,  there  are  both  in-flight  and  laboratory  procedures  and  experiments 
that  can  be  utilized  to  investigate  this  hypothesis  and  its  potential  impact  further; 
for  example,  sampling  simulation  experiments  through  a  flow  tube  in  the  laboratory; 
in-flight  operational  changes  such  as  intermittent  thermal  desorption  of  wall  deposits 
and  re-equilibration  prior  to  sampling,  and/or  sampling  at  lower  altitudes  first  lor 
NO^;  and  so  on.  However,  the  feasibility  of  carrying  out  any  of  these  efforts  must 
be  evaluated  within  the  context  of  current  and  projected  activity  of  the  sampling 
program. 

Jr  addition  to  the  two  preceding  scenarios  involving  the  possible  effects  of 
heterogeneous  phenomena  at  the  sampling  tube  walls  on  the  total  cryogenic  sample 
co.  lposition,  there  are  probably  others  that  can  be  formulated  as  well  that  are  in 
qualitative  agreement  with  the  known  data.  Without  additional  experimental  evidence 
regarding  the  precise  nature  of  the  processes  at  the  tube  wall,  however,  these  must 
remain  conjectures  only. 

4.3  Effect  of  Tube  Wall  Temperature 

In  considering  factors  that  could  significantly  change  the  magnitude  of  the  rate 
constant  values  assumed,  the  possibility  of  radiative  solar  heating  if  the  Men-iwav 

sample  tubing  was  examined.  Although  no  temperature  measurements  ot  the 

* 

sampling  tubing  were  made  directly,  both  the  air  temperature  ana  the  external 
surface  temperature  of  the  TKIW'AS  unit  (at  a  point  midway  between  the  top  an  i 
bottom)  were  recorded.  This  surface  temperature  should  be  roughiv  equivalent  to 
the  sampling  tube  wall  temperature.  As  shown  in  Table  a,  the  temperature  lifter  • 
cnees  between  the  surface  and  air  were  both  positive  <«I°  C  maximum  >  and 
negative  (-f>’t',  maximum).  These  differences  in  temperature  should  not  lie  large 
enough  to  cause  anv  serious  increase  (or  decrease  I  in  air  sample  temperature  due 
to  heat  transfer  through  the  wall  of  the  sampling  tubing.  An  estimate  of  the  potential 
magnitude  of  this  effect  follows. 

Assuming  a  constant  tube  wall  temperature,  T  ,  an  energy  balance  on  the 
I  sampled  air  is; 

2rr  HI.  h  AT  W  C  IT.,  -  T.i  1 3  1  i 

)  avg  p  2  1 

t 

\  _ 

* 

Note  added  in  proof;  Since  this  report  was  written,  a  balloon  flight  oackagc  using 
the  same  sampling  tube  was  fitted  with  thermistors  placed  at  various  locations  and 
flown  in  the  lower  stratosphere.  The  maximum  surface  temperature  throughout 
.  the  13-13  km  sampling  range  was  -17’C.  Thus,  dearie  all  the  conclusions 

arrived  at  in  this  report  are  supported,  and,  in  Ice  i,  t-  av  actually  he  conserval  i\ 
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laminar  Hon  wall  heat  1  ranslV  r  eoefficient,  eal  /  chi'"*  see*  K , 
loi.*;*  r  i  1  hin  i  >  mean  l«  m  >.  t »« •  rature  -  ii  t'f'erenre  between  the  mi  be*  wall  and 
average  hulk  ten -pernture  at  the  inlet  anil  the  sampling  point, 
air  mass  Ihnv  rale,  ^f.scc, 
speeific  heat  <  apa<  itv  >f  the  air,  cal  t»*K, 
hulk  temperature1  at  sampling  point., 
hulk  temperature  at  inlet. 


The  \usselt  number  (2hK/k  ;  k  =  thermal  conductivity  of  air)  for  laminar  tube  flow 
vv  ith  constant  wall  temperature  is  well  known  (for  example,  see  J$in!  et  al°,  p.  406). 
l  or  the  three  sampling  altituiles  considered  here  and  the  How  conditions  in  Sta  ¬ 
tion  2.  1 : 
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where  I'  T,  Dl.h'tt 
presented  in  Table 


Altitude 


l*  ■  Is  inn  the  extreme  values  of  the  temperature  differences 
P 

a,  together  with  the  basic  data  from  Section  2.  1: 
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Another  method  of  assessing  the  magnitude  of  t  .e  heat  transfer  effect  is  to 
estimate  the  tube  wall  temperature  required  to  increase  the  air  sample  ten  nera- 
ture  bv  10K; 


Altitude  Wall  Temperature  Required 

(km)  to  Heat  Air  Sample  10k  AT  Above  Ambient 

(K  )  ( K i 


la 

72  1.6  (446,  5“  Cl 

505. 

20 

40  2.  (218.  H 0  (  ) 

275.  4 

TO 

280.  (7,T°C> 

40.  2 

T  hus,  even  at  TO  km,  there  is  no  reason  to  suspect  an  air  sample  temperature 
increase  (or  decrease!  even  approaching  IOK;  that  is,  the  4><.  2K  temperature  differ¬ 
ence  required  is  significantly  greater  than  the  KK  or  so  suspected.  The  conclu¬ 
sion  of  this  analysis  is  that  the  effect  if  temperature  on  processes  within  the  sampling 
ubc  is  completely  negligible.  In  addition,  since  the  tube  wall  temperatures  are 
expected  to  be  within  IOK  of  ambient,  heterogeneous  interactions  should  also  occur 
at  approximately  ambient  temperatures. 


5.  SUMMARY  ANI)  CONCLUSIONS 

A  numerical  model  of  a  sampling  tube  used  to  conduct  stratospheric  air  to  the 
sampling  point  of  a  balloon-borne  cryogenic  whole  air  sampler,  has  been  developed 


♦ 

The  effect  of  altitude  is  reflected  bv  the  ratio,  p.  Both  the  heat  transfer  coefficient, 
h,  and  the  mass  flow  in  the  sampling  tube,  W  ,  decrease  with  altitude;  however,  h 
decreases  more  slowly  than  W. 


i 

I 
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T  6 


in  order  to  assess  the  potential  effects  of  passage  through  the  tube  on  alteration  of 
species  mixing  ratios  from  those  in  the  ambient.  This  model  is  based  on  the  applica¬ 
tion  of  a  general  purpose,  homogeneous,  multireaction  chemical  kinetics  code 
(CHKMSMN)  to  a  3  1  reaction  (4!.*  reaction  with  heterogeneous  wall  loss),  23  species 
(42  species  with  adsorbed  or  "lost"  constituents)  stratospheric  kinetic  model, 
modified  to  take  into  account  heterogeneous  interactions  with  the  tube  wall.  Ilssen- 
lially,  the  heterogeneous  loss  rates  due  to  radial  diffusion  were  transformed  into 
equivalent  pseudo -homogeneous  reactions.  This  technique  exploits  the  difference 
in  time  scales  for  radial  diffusion  and  convection  in  the  sampling  tube.  The  re¬ 
sultant  model  represents  a  useful  general  tool  lor  estimating  the  potential  absolute 
and  relative  effects  of  the  tube  wall  and  homogeneous  chemical  reactions  in  altering 
the  composition  of  complex  flowing  and  reacting  mixtures  in  atmospheric  or  labora¬ 
tory  applications. 

The  effects  of  the  homogeneous  gas  phase  rhemistrv  and  heterogeneous  wall 
loss  were  estimated  at  three  different  altitudes:  15,  20,  and  20  km.  It  was  found 
that  the  gas  phase  chemistry  was  generally  incapable  of  appreciably  altering  the 
sampled  composition  from  that  in  the  ambient.  This  result  is  primarily  due  to  the 
short  average  residence  time  in  the  sampling  tube.  On  the  other  hand,  depending 
>n  the  particular  set  of  assumptions  used,  heterogeneous  wall  loss  was  found 
capable  of  imposing  significant  alterations  on  the  sampled  gas  composition,  parti¬ 
cularly  at  altitudes  above  about  20  km.  due  to  fast  radial  diffusion.  The  tube  wall 
temperature  was  shown  to  have  a  completely  negligible  effect  on  the  homogeneous 
gas  phase  kinetics. 

Currently,  definitive  quantitative  assessments  of  die  impact  of  the  heterogeneous 
interactions  at  the  tube  wall  on  the  final  sampled  coyp position  are  hampered  bv  lack 
if  knowledge  concerning  the  precise  nature  of  the  complex  processes  (of  all  those 
possible)  that  actually  occur  to  a  significant  extent  at  stratospheric  conditions  and 
compositions.  In  this  regard,  the  behavior  of  the  tube  wall  with  respect  to  anv 
single  species  may  range  from  that  of  a  simple  sink  (of  varying  effectiveness,  as 
reflected  bv  the  capture  coefficient),  to  that  of  a  complex  chromatographic  chemical 
reactor.  However,  two  scenarios  involving  heterogeneous  wall  effects  were  sug¬ 
gested  that  at  least  qualitatively  explain  some  trends  observed  in  comparisons  of 
v  arious  other  data  and  model  predictions  w  ith  the  Al  t'.  I.  ervosampler  data.  The 
"oiltaassing’  hypothesis  in  particular,  whereby  NO  species  collected  at  high  alti¬ 
tudes  could  he  released  at  lower  altitudes,  seems  deserving  of  further  study. 

In  order  to  improve  quantitative  predictive  techniques  for  the  effects  of  sampling 
network  walls  on  sampled  compositions,  good  experimental  data  are  needed  to  he 
used  in  conjunction  with  a  practical  numerical  model.  This  report  provides  an 
approach  and  the  rudiments  of  an  appropriate  model.  However,  t Ho  c  ontrolling 
heterogeneous  physicochemical  processes  for  species  of  interest  at  st ratnspheri r 
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Appendix  A 

Newton  s  Method  Solution  for  [Eq.  (23)] 


Application  of  Newton's  method  to  find  the  roots  of  Eq.  (23)  can  be  expressed 


<i+  •  >  _  „  <»>  .  ‘‘N/o'.a  (i,> 


(AD 


-  x  =  a  j"  -  g(«  t  '/b  t 

where  (i)  represents  the  current  estimate  fora  (,  gta  is  the  value  of  the  function, 
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wherein  the  arguments  of  all  the  Bessel  functions  are  (a  ,U  D  !<>• 

Iterative  application  of  Eq.  (Al)  to  convergence,  using  some  predetermined 
criterion,  will  yield  the  value  of  a  j  for  the  parameters  f,  d.  and  It. 
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HPRiME  h  /  7  4 
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Subroutine  KPRIME 
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